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ABSTRACT 

The young X-ray and gamma-ray-bright supernova remnant RX J1713.7— 3946 
(SNR G347.3— 0.5) is believed to be associated with molecular cores that lie within re- 
gions of the most intense TeV emission. Using the Mopra telescope, four of the densest 
cores were observed using high-critical density tracers such as CS(J=l-0, J=2-l) and its 
isotopologue counterparts, NH3(1,1) and (2,2) inversion transitions and N2H+(J=l-0) 
emission, confirming the presence of dense gas >10*cm~^ in the region. The mass 
estimates for CoreC range from 40Mq (from CS(J=l-0)) to SOM© (from NH3 and 
N2H+), an order of magnitude smaller than published mass estimates from CO(J=1-0) 
observations. 

We also modelled the energy-dependent diffusion of cosmic-ray protons acceler- 
ated by RX J1713.7— 3946 into CoreC, approximating the core with average density 
and magnetic field values. We find that for considerably suppressed diffusion coeffi- 
cients (factors X =10"'^ down to 10^^ the galactic average), low energy cosmic-rays 
can be prevented from entering the inner core region. Such an effect could lead to 
characteristic spectral behaviour in the GeV to TeV gamma-ray and niulti-keV X-ray 
fluxes across the core. These features may be measurable with future gamma-ray and 
multi-keV telescopes offering arcminute or better angular resolution, and can be a 
novel way to understand the level of cosmic-ray acceleration in RX J1713.7— 3946 and 
the transport properties of cosmic-rays in the dense molecular cores. 

Key words: diffusion - molecular data - supernovae: individual: RX J1713.7— 3946 - 
ISM: clouds - cosmic rays - gamma rays: ISM. 



1 INTRODUCTION 

The origin of galactic cosmic rays (CRs) is mysterious, 
but popular theories attribute observed fluxes to first or- 
der Fermi acceler ation in the shocks of sup ernova rem- 
nants (SNRs) fee. iBl andford fc Ostrikeil (HHl)). One such 
potential CR-a ccelerator is RX J 1713.7— 3946, a young, 
~1600 year-old l|Wang et al~lll997l) SNR that exhibits shell- 
like properties at both X-ray llPfeffermann et al. I [l996l: 
Kovarna et al. 19971: iLazendid 2003 : Cassam-Chenai et al. I 



2004 iTanaka et al. I 12008^: 'Acero et al. ' '2009*1 and TeV 
gamma-ray energies lAharoni an et al. 2006 , 2007). 

A void in molecular gas, most likely created by 
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the RX J1713.7— 3946 p rogenitor star-winds (jlnoue et al. I 
l2012l : iFukui et al.1l20l2l ). has previously been observed in 
CO(J=1-0) an d CO(J=2-l) with the Nanten and Nanten2 
4 m t elescope (|Fukui et al. I l2003l : iMoriguchi et al. I l2005l : 
lFukui2008 l. and more rece ntly, HI in the Southern Galac- 
tic Plane Survey (SGPS) (|McClure-Grifiiths et al. I l2005l : 
iFukui et al. II2012I '). Three dense clumps. Cores A, B and C 
(see Figure[l]), to the West of this void, and a fourth clump, 
CoreD, in the remnant's northern boundary are coincident 
with regions of high TeV flux. The beam FWHM of ~10' for 
the detection of gamma-rays with HESS would not be able 
to resolve TeV features of the scale of these molecular cores 
(-!')• 

Three of the four stated cores (Cores A, C and D) 
feature coincident infrared sources commonly associated 
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with star formation l|Moriguchi et al. 1 12005| ). Core C dis- 
plays a bipolar nature in CO(J=3-2) and CO(J=4-3) that 
may be attributed to protostellar act ivity or outflows 
l|Moriguchi et al. II2005I : ISano et al. II2OI0I ). 

RX J1713.7— 3946 has a well-defined high energy struc- 
ture that is shell-like at both keV X-ray and TeV gamma-ray 
energies (|Aharonian et al. I [iooil '). The TeV emission con- 
tours displayed in Figure [TJi reveal a clear shell-like struc- 
ture. The highest fluxes originate from an arc-shaped region 
in the northern, north-western and w estern parts of the rem- 
nant . The keV XMM-Newton image (jCassam-Chenai et al. I 
l2004h consists of extended emission regions that form what 
appears to be two distinct arced shock fronts in the north, 
north-west and western parts of Figure [1]d. Two separate 
shock structures seen at keV energies towards the north- 
west, believed to be dominated by synchrotron radiation 
from electrons within the RX J1713.7— 3946 shock, have 
been considered as an outwards- moving and a re flected, 
inwards- moving shock (|Cassam-Chenai et al. |[2004h . Vari- 
ous flux enhancements and filaments are evident along these 
fronts. 

The RXJ1713. 7-3946 keV emission was first thought 
to originate from a distance of ~6kpc due to an assumed 
association with molecular gas at line-of-sight velo cities, 

vlsr ' 70kms"^ to -90kms"^ ijSlane et ainil999l), but 

later CO observations and discussion bv'Fuku i et al. I ()2003l ) 
placed the remnant inside a void in molecular gas at ~1 kpc 
distance (vlsr ~ — 10kms~^). This distance was consistent 
with an initial estimate ba sed on an X-ray spectral fit and 
a Sedov model llfeffer ma nn et al. Ill9 96). which was itseli' 



consistent with the SNR age (jWanget al. Ill997h . 

Interestingly, keV X-ray emission peaks (Figure [TJj) 
seem to show some degree of anticorrelation with molecu- 
lar gas peaks and cores at vlsr ~ — lOkms"^, as traced 
by CO (Figure [T}d). Two synchrotron intensity peaks on the 
border of the densest central region of Core C are sugges- 
tive o f compressions associated with the shock (jSano et al. I 
I2OIOI ). while a dip in non-thermal emission from CoreC's 
centre may imply that the electron population is unable 
to diffuse into the core and/or that keV photons can- 
not escape. A similar scenario may also be applicable to 
CoreD, with an X-ray intensity peak also coincident with 
CoreD's boundary. Additional gas-keV association may be 
seen in a line running along an X-ray filament on the 
eastern border of Core A, Point ABm (Midpoint of Cores 
A and B) and CoreB. These structures possibly suggest 
gas swept -up by the RX J17 13.7— 3946 shock or progeni- 
tor wind. Ilnoue et al. I \2Q1± modelled a SN shock-cloud 
interaction and outlined the resultant high energy emis- 
sion characteristics that were common to both their model 
and RX J1713.7— 3946. These included parsec-scale keV X- 
ray correlation with CO emission, subparsec X-ray anti- 
correlation with CO emission peaks, short-term variabil- 
ity of bright X-ray emission, and stronger TeV gamma- 
ray emission in regions with stronger CO emission. De- 
spite this strong evidence fo r an association, a survey by 
ICassam-Chenai et al. I l|2004h that t argeted shock-tracing 
1720 MHz OH maser emission ijWardle 1999) recorded no de- 
tections towards this gas complex. This lack of an indication 
of a molecular chemistry consistent with a shock may point 
towards a scenario where the gas does not lie within the 
SNR shell, but given the weight of evidence to the contrary 



(outlined in this se ction), the most probable ex planations 
are, as suggested bv lCassam- Chenai et al"l l|2004h , that the 
density is not sufficiently large to cause OH maser emission 
or that a non-dissociative C-type shock, like that required, 

does not exist in this region. 

Fermi-LAT observations (|Abdo et al. I l201lf ) have 
recently shown that RX J1713. 7-3946 exhibits a low 
but hard-spectrum flux of 1-10 GeV gamma-ray emis- 
sion, not predicted by early hadronic emission models 
(ie. pion-producing proton- proton collisions), but con- 
sistent with previously- published lepton-dominated 
gamma-ray models (ie. inverse Compto n scattering 
of T e V electrons) of RXJ1 71 3. 7-3946 JPorter et al. I 



20061: lAharonian et all l2007l: iBerezhko fc Vollj bOloT 
Ellison et al. II2OI0I : IZirakashvih fc AharonianI I2OI0I I. How 



ever, a hadronic component is possible towards dense 
cloud clumps ( Zirakashvili & Aharonian 2010), and per- 
haps more globally if one considers an inhomogeneous 
ISM (jlnoue et al. 1120121 ) into which the SNR shock has 
expanded. Further support for a global hadronic component 
may come from consi deration of the mol ecular and atomic 
ISM gas together (jFukui et al. I I2OI2I ). Another novel 
way to discern the leptonic and/or hadronic nature of 
the gamma-ray emission is to make use of the potential 
energy- depend ent diffusion o f CRs into dense cloud clumps 
or cores (e.g. ICabici et al.1 l|2007 ')). which may lead to 
characteristic features in the spectrum and spatial distri- 
bution of secondary gamma-rays and synchrotron X-rays 
from secondary electrons. In §[6] we therefore looked at a 
simplistic model of the energy-dependent diffusion of CRs 
into a dense molecular core similar to Core C, in an effort to 
characterise the penetration depth of CRs in the core and 
qualitatively discuss the implications for the gamma-ray 
and secondary X-ray synchrotron emission. 

CO observations, such as those already obtained to- 
wards RX J1713.7— 3946, are ideal for tracing moderately 
dense gas (~10'^ cm~'^), but they do not necessarily probe 
well dense regions and cores (10^ cm~"^) that may play an 
important role in the transport and interactions of high en- 
ergy particles. Dense gas may also influence SNR shock 
propagation. To gain a more complete picture, we per- 
formed 7 mm observations of high-critical density molecules 
to probe the inner regions of cores and search for evidence 
of ISM disruption by the SNR shocks and star-formation. 
Towards this goal, we used the Mopra 22 m single dish 
radiotelescope in a 7 mm survey to observe bands con- 
taining the d ense gas trace r CS(J = l-0), the shock tracer 
SiO( J =l-0') JSchilke et al. I 1 19971 : iMartin-Pintado et al. I 
I2OO0I : ICusdorf et al. ' '2008a'^ ) and sta r- formation tracers 
CH3OH and HC3N (van Dishoek fc Blake 199a ). We chose 
the western region of RX J1713.7— 3946 due to the evidence 
of shock-interactions with the molecular gas. 

Additionally, a 12 mm observation of NH3 inversion 
transitions and archival 3mm observations of CS(J=2-1) 
provided further information towards Core C. 



2 OBSERVATIONS 

An 18'xl8' region (see Figure[l| centred on (Z,6) = (346.991, 
-0.408) that encompasses Cores A, B and C, was mapped 
by the 22 m Mopra telescope in the 7 mm waveband on the 
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Figure 1. a)(lcft) Nanten2 CO{J=2-l) inte grated-intensity image (vlsr = -18 to Okms-l) l|Sano et al. 1120101 ^ with HESS >1.4TeV 
gamma-ray excess count contours overlaid fAharonian et al. 2006). Contour units are excess counts from 20 to 50 in intervals of 5. 
All core names except ABm were assigned by Moriguchi ct al. (2005). The dashed 18'xl8' box marks the area scanned in 7mm 
wavelengths by Mopra (Figure |2]l. b)(right) XMM-Newton 0.2-12 keV X-ray image with Nanten2 CO(J=2-l) {vlsr=-18 to Okms'^) 
contours. CO(J=2-l) contours span 5 K km s~^ to 40 K km s~^ in increments of 5 K km s~^. The XMM-Newton image has been exposure- 
corrected and smoothed with a Gaussian of FWHM=30". Two compact objects are indicated: IWGA J1 713.4— 3949, which has been 
suggested to be associated with R X J1713.7— 3946 and IWGA J1714.4— 3945, which is likely unassociated jSlane et al. Ill999l : lLazendi3 
I2OO3I: ICassam-Chenai et al. II2004I) . 



Table 1. The window set-up for the Mopra Spectrometer, MOPS, 
at 7 mm. The centre frequency, targeted molecular line, targeted 



frequency and total mapping noise (Trms 


) are displayed. 


Centre 


Molecular 


Line 


Map 


Frequency 


Emission Line 


Frequency 


TrmS 


(GHz) 




(GHz) 


(K/ch) 


42.310 


30SiO(J=l-0,v=0) 


42.373365 


~0.07 


42.500 


SiO(J=l-0,v=3) 


42.519373 


~0.07 


42.840 


SiO(J=l-0,v=2) 


42.820582 


~0.07 




29SiO(J=l-0,v=0) 


42.879922 




43.125 


SiO(J=l-0,v=l) 


43.122079 


~0.08 


43.395 


SiO(J=l-0,v=0) 


43.423864 


~0.08 


44.085 


CH3OH(7(0)-6(l) A-H-) 


44.069476 


-0.08 


45.125 


HC7N(J=40-39) 


45.119064 


~0.09 


45.255 


HC5N(J=17-16) 


45.26475 


~0.09 


45.465 


HC3N(J=5-4,F=5-4) 


45.488839 


-0.09 


46.225 


i='CS(J=l-0) 


46.24758 


-0.09 


47.945 


HC5N(J=16-15) 


47.927275 


-0.12 


48.225 


C34S(J=l-0) 


48.206946 


-0.12 


48.635 


OCS(J=4-3) 


48.651604 


-0.13 


48.975 


CS(J=l-0) 


48.990957 


-0.12 



nights of the 22nd of April 2009, 23rd of April 2009, 11th 
of April 2010 and the 21st of April 2010. The Mopra spec- 
trometer, MOPS, was employed and is capable of recording 
sixteen tunable, 4096-channel (137.5MHz) bands simultane- 
ously when in 'zoom' mode, as used here. A list of mea- 
sured frequency bands, targeted molecular transitions and 
achieved Trms levels are shown in Table [1] 

In addition to mapping, 7 mm deep ON-OFF switched 



pointings (pointed observations) were taken on selected re- 
gions on the nights of the 23rd of April 2009, 16th of March 
2010, 17th of March 2010 and 18th of March 2010, and a 
12mm pointed observation was taken on the night of the 
10th of April 2010. Archival 3mm Mopra data on Core C 
were utilised in the analysis and were taken on the 30th of 
April, 2007 in the MOPS wide-band mode (8.3 GHz band). 

The H2Q Southern Galactic Plane Survey (HOPS) 
(|Walsh et al. I [2008) surveyed the entire Galactic Longitu- 
dinal extent of RX J1713.7— 3946 down to a Galactic Lat- 
itude of -0.5°, reaching a noise level of Trms ~0. 25 K/ch. 
From these data, only a weak NH3(1,1) detection towards 
Core C was found. We took deeper 12mm pointed observa- 
tions towards this region in response to the HOPS detec- 
tion, revealing NH3(1,1) and NH3(2,2) emission of a suffi- 
cient signal-noise ratio to extract some gas parameters. 

The velocity resolutions of the 7 and 12 mm zoom- 
mode and 3 mm wideband-mode data are — 0.2kms~^, 
~0.4kms~^ and ~lkms~^ respectively. The beam FWHM 
of Mopra at 3, 7 and 12mm are 36±3", 59.4±2.4" and 
123±18", respectively, and the pointing accuracies are ~6". 
The achieved Trms values for 7 mm maps and individual 3, 
7 and 12 mm pointings are stated in Tables [1] and (2] 

On-The-Fly (OTF) mapping and pointed observation 
data were reduced and analysed using the ATNF anal- 
ysis programs, Livedata, Gridzilla, Kvis, Miriad and 

asafQ. 

Livedata was used to calibrate each OTF-map scan 



^ See |http://www.atnf.csiro.au/computing/software/| 
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(row/column) using the intermittently measured back- 
ground as a reference. A linear baseline-subtraction was 
also applied. Gridzilla then combined corresponding fre- 
quency bands of multiple OTF-mapping runs into 16 three- 
dimensional data cubes, converting frequencies into line-of- 
sight velocities. Data were weighted by the Mopra system 
temperature and smoothed in the Galactic I — b plane using 
a Gaussian of FWHM 1.25'. 

MiRlAD was used to correct for the efficiency of the 
instrument l|Urquhart et al. 1|2010| ) for map data and create 
line-of-sight velocity-integrated intensity images (moment 0) 
from data cubes. 

ASAP was used to analyse pointed observation data. 
Data were time-averaged, weighted by the system temper- 
ature and had fitted polynomial baselines subtracted. Like 
mapping data, deep pointing spectra were corrected for th e 
Mopra efficiency (|Ladd et al. ..2005. : .Urguhart et al.~ll2010l V 



3 LINE DETECTIONS 

This investigation involved 6 different species of molecule 
and all detections are displayed in Table [2] Our 7 mm sur- 
vey found the dense gas tracers CS, C''''S and ^''CS in 
the J=l-0 transition and HC3N(J=5-4) (C^''S(J=l-0) and 
HC3N(J=4-5) emission maps are in §|Al]|. Various detec- 
tions of CH3OH are indicative that these cores have a warm 
chemistry. Deep 12 mm observations measured two lines of 
rotational NH3 emission. In addition to these molecules, 
archival 3 mm data revealed detections of the J=2-l tran- 
sitions of CS and C'^'*S, transitions of SO and N2H''' and 
further detections of CH3OH. 



3.1 CS(J=l-0) Detections 

Figure [5] is a CS(J=l-0) map of the CoreA-B-C region of 
RXJ1713. 7-3946. CS(J=l-0) emission can be seen to be 
originating from the two most C0(J=2-l)-intense regions. 
Cores A and C at line-of-sight-velocities of ~ — lOkms"^ 
and ~ — 12kms~^, respectively, consistent with CO emis- 
sion line-of-sight-velocities. The CoreC CS(J=l-0) emission 
is coincident with the Core C infrared source, but the Core A 
CS(J=l-0) emission peak is offset from the Core A infrared 
source. This discrepancy of ~l', similar to the size of the 
7mm beam FWHM, is not evident for CO(J=2-l) emission. 

Figures |3 H O [3 H |8] and H show CS(J=l-0) spec- 
tra from pointed observations towards Cores A, B, C 
and D, and Point ABm . Corresponding CO(J=1-0) spectra 
jMoriguchi et al. II2005I ) are also shown in these figures, but 
the CO beam FWHM of ~180", corresponding to ~9x the 
7 mm solid angle, means that neighbouring gas is probably 
included in the beam-average. 

Core A exhibits symmetric narrow line CS(J=l-0) emis- 
sion, of Avfwhm ~1.25kms~^, suggesting that the in- 
ner, dense core is not under the infiuence of an exterior 
shock. Figure (TJa shows that the forward shock is coin- 
cident with the outskirts of Core A such that the inner 
CS(J=l-0)-emitting region probably remains unaffected by 
the RX J1713.7— 3946 shock, consistent with our observa- 
tions. 

CoreC displays the strongest CS(J=l-0) emission of 
this survey, implying that this region probably harbours 




Galactic Longitude (cleg) 

Figure 2. Mopra CS(J=l-0) integrated-intensity image (vlsr=- 
13.5 to -9.5kms-i) with Nanten2 CO{J=2-l) (vlsr=-18 to 
Okms"'^) contours (black) and two 8 /xm point sources indicated 
by crosses. CO(J=2-l) contours span 5Kkms~^ to 40Kkms~^ 
in increments of 5Kkms~^. CS(J=l-0) contours (white) are 
also shown, spanning from 0.4 to 2.0Kkms~^ in increments of 
0.4Kkms-l. The beam FWHM of Mopra at 49GHz is displayed 
in the bottom left corner. 
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Figure 3. Core A: CS(J=l-0) spectrum (top) and CO(J=1-0) 
spectrum (bottom) towards Core A. The position of Core A is that 
displayed in Table[2] 
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Figure 4. CoreB: CS(J=l-0) spectrum (top) and CO(J=1-0) 
spectrum (bottom) towards Core B. The position of Core B is that 
displayed in Table|2] 
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Table 2. Detected molecular transitions from pointed observations. Velocity of peak, vlsr, peak intensity, Tpoak, and FWHM, Aufwhm, 
were found by fitting Gaussians before de convolving with the MOPS velocity resolution. Displayed values include the beam efficiencies 
llUrguhart et al. II2OI0I : iLadd et al.1l2005l) . after a linear baseline subtraction. Statistical uncertain ties are show n, whereas systematic 
uncer tainties are ~7%, ~2.5% and ~5% for the 3mm, 7mm and 12 mm calibration, respectively llLadd et al. I 2005; Urquha rt et al. I 
l201(tl . Band noise, Trms, integrated intensity, jT^jdv, possible counterparts and 12 and 100 /im IRAS flux, F12 and Fioo, respectively, 
(where applicable) are also displayed. 
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(K/ch) 


(km s""'") 


(K) 


(kms-l) 


(Kkms-1) 


fFio/Finn (Jv)] 


Core A 


CS(J=l-0) 


0.06 


-9.82±0.02 


0.92±0.04 


1.25±0.06 


17.0±0.9 


IRAS 17082-3955 


(346.93°, -0.30°) 


HC3N(J=5-4,F=5-4) 


0.04 


-9.76±0.04 


0.29±0.02 


1.07±0.11 


4.6±0.6 


[5.4/138] 




SiO(J=l-0) 


0.03 


_ 


_ 


_ 






CoreB 


CS(J=l-0) 


0.06 


- 


- 


- 


- 




(346.93°, -0.50°) 


SiO(J=l-0) 


0.03 


- 


- 


- 


- 
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0.10 


J- -L . U^^U. UO 






i\A ^-1-9 

U'i. O^Z . U 


TRA^ 17nSQ— 'iQ'^l 


(347.08°, -0.40°) 


C34s(J=l-0) 


0.04 


-11.83±0.04 


0.37±0.02 


1.40±0.09 


7.7±0.7 


[4.4/234] 




C34S(J=2-1) 


0.07 


-11.78±0.20 


0.34±0.06 


1.89±0.49 


9.5 ± 1.4 






i3cS(J=l-0) 


0.03 


-11.70±0.10 


0.10±0.02 


0.88±0.22 


1.3 ± 0.5 






HC3N(J=5-4,F=5-4) 


0.02 


-11.65±0.03 


0.32±0.01 


1.74±0.08 


8.2 ± 0.6 
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CH30H(2(-1)-1(-1)E) 


0.07 


-11.75±0.17 


0.33±0.13 


0.99±0.90 


4.8 ± 1.7 






CH3OH(2(0)-l(0)A++) 




-11.86±0.14 


0.43±0.08 


1.44±0.46 


9.2 ± 1.5 






SO(2,3-l,2) 


0.07 


-11.77±0.09 


0.65±0.07 


1.74±0.25 


16.7 ± 1.4 






N2H+(J=l-0) Fi=2-1" 


0.08 


-11.61±0.03 


1.14±0.03 


2.18±0.08 


36.8 ± 1.1 






Fi=0-1 


// 


// 


0.40±0.04 


1.58±0.23 


9.4 ± 1.1 






Fi = l-1 


// 


// 


0.76±0.03 


0.87±0.04 


9.8 ± 0.6 






NH3((1,1)-(1,1)) F-^F" 


0.03 


-12.04±0.04 


0.34±0.02 


1.46±0.09 


7.3 ± 0.7 






F-s>F-l'' 


n 


// 


0.07±0.02 


// 


1.5 ± 0.5 








n 


// 


0.13±0.02 


// 


2.8 ± 0.5 






F-5>F-fl= 


ti 


// 


0.13±0.02 


// 


2.8 ± 0.5 








ti 


// 


0.12±0.02 


// 


2.6 ± 0.5 






NH3((2,2)-(2,2)) F^F" 


II 


-11.69±0.10 


0.16±0.02 


2.02±0.26 


4.8 ± 0.8 






SiO(J=l-0) 


0.03 












CoreD 


CS(J=l-0) 


0.10 


-9.06±0.11 


0.47±0.04 


2.50±0.29 


17.4 ± 1.3 


IRAS 17078-3927 


(347.30°, 0.00°) 




0.10 


-71.25±0.46 


0.15±0.04 


3.38±0.36 


7.5 ± 1.2 


[2.0/739] 




SiO(J=l-0) 


0.06 












Point ABm 


CS(J=l-0) 


0.04 


-8.89±0.43 


0.09±0.02 


4 ±1 


5.3 ± 1.1 




(346. 93°, -0.38°) 


SiO(J=l-0) 


0.06 













"■Centre line, ^ Outer satellite lines, Inner satellites lines, 



the larsest mass dense core. CO(J=4-3) observations by 
ISano et al. I (|2010l ) suggested that Core C contains a bipolar 
outflow, but we find no obvious asymmetries in the CS(J=1- 
0) emission line profile or map data. We do, however find a 
significant difference between the red and blue-shifted sides 
of CS(J=2-1) emission. For deep switched pointing spectra, 
fitted gaussians were subtracted from the observed spectra 
of CS(J=1-0,J=2-1) emission, and the resultant spectra were 
integrated over velocity ranges chosen, by eye, to represent 
red and blue-shifted line wings. Integrating CS(J=2-1) emis- 
sion over 4.1 kms~^-wide bands 3.0kms~^ either side of the 
peak emission, highlights a 2-2. 5(j intensity difference be- 
tween red and blue-shifted sides, possibly trac ing manifes- 
tation s of the Core C bipolar outflow seen by ISano et al. I 
in CO emission lines. No significant asymetry (<lcr) 
is viewed in CS(J=l-0) emission spectra, while in the plane 
of the sky, there is no significant offset between the red and 
blue-shifted CS(J=l-0) emission peaks. 

In the northern part of the remnant, CoreD exhibits a 
CS(J=l-0) emission line with a FWHM of 2.5±0.3kms-\ 



slightly larger than that for Cores A and C. This might 
be expected given that Core D harbours an IRAS point 
source more luminous than t hat of the other cores asso- 
ciated with RXJ1713. 7-3946 (|Moriguchi et al. II2005I ). but 
CoreD also lies on the northern shock front, another possi- 
ble source of disruption. Towards CoreD, a second peak at 
Vlsr. =~ -70kms-i (~6kpc) was detected. This Norma- 
arm gas, initially suggested bv lSlane et a Tl J 19991) to be as- 
sociated with RX J1713.7— 3946, is probably unassociated, 
given the preference for a Ikpc distance (see iJT|. 

Deep observations on a location at the midpoint of 
Cores A and B, Point ABm, revealed weak and moderately 
broad CS(J=l-0) emission of FWHM, Ai;fwhm=4±1, pos- 
sibly indicating the existence of a passing shock, as there 
is no evidence to support star-formation towards this lo- 
cation. This scenario is qualitatively supported by the po- 
sitional coincidence of Point ABm with the western shock- 
front (see but the low signal to noise ratio means that 
any conclusions drawn from this are low-confidence. Alter- 
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Figure 5. Core C: CS isotopologuc spectra towards CoreC. The 
position of Core C is that displayed in Tabled 
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Figure 6. C ore C : Various CO sp e ctra to wards Core C 
jPukui et al. I |2003|: iMoriguchi et al. I |2005|: iFukuil |2008|; 
ISano et al. I |2010|) . The position of CoreC is that displayed 
in TableO 



natively, small-scale unresolved clumps (perhaps from star- 
formation) may be responsible for the detection. 



3.2 Other Molecular species 

No emission from the shock tracer, SiO, was seen in maps 
or pointed observations. The RX J1713.7— 3946 shock-front 
likely exceeds 1000 km s~^, which may dissociate molecules 
in the region, but also initiate no n-dissociative sub s idiary 
shocks in dense e discussed bv lUchivama et al. I (l20ld ) 
for older SNRs. Given the existence of dense gas in the re- 
gion, as indicated by this study, one might indeed expect 
such secondary shocks if the gas lies within the SNR-shell. 
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Figure 7. Core C: Various spectra towards Core C. The position 
of Core C is that displayed in Table[2] 
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Figure 8. CoreD: CS(J=l-0) spectrum (top) and CO(J=1-0) 
spectrum (bottom) towards CoreD. CS(J=l-0) emission has been 
box-car-smoothed over 6 channels. The position of Core D is that 
displayed in Tablc[2] 



In such scenarios, it may be expected that Si is sputtered 
from dust grains, increasing the gas-phase abundance of 
SiO in the post-shock region (see for example ISchilke et al. I 
(|1997M . As we do not detect SiO towards this region of 
RXJ1713. 7-3946, we conclude that, if the SN shock is in- 
deed interacting with the surveyed gas as indicated by high 
energy emissions (see §[l|, it has not produced SiO at de- 
tectable levels. 
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Figure 9. Point ABm: CS{J=l-0) spectrum (top) and CO(J=l- 
0) spectrum (bottom) towards Point ABm. CS(J=l-0) emission 
has been box-car-smoothed over 6 channels. The position of 
Point ABm is that displayed in Table|2] 



Detections of the organic molecules HC3N and CH3OH 
were recorded from CoreC and HC3N from Core A. CH3OH 
is believed to be the product of grain-surface reactions. Weak 
emission from this molecule suggests that conditions within 
the core are capable of evaporating a significant amount of 
the ice mantles of grains into gas phase. 

N2H"'" and NH3 were both detected inside CoreC. NH3 
is generally thought to be a tracer of cool, quiescent gas, 
as its main reaction path ways do not require high temper- 
atures (|Le BourlotI Il99l} ). The main reaction pathway of 
N2H"'" is thro ugh molecular nitrogen interacting with 
iHotzel et al. . 2004), so N2H''" can be thought as a tracer of 
molecular nitrogen. 



4 MILLIMETRE EMISSION LINE ANALYSES 

Local Thermodynamic Equilibrium (LTE) analyses us- 
ing NH3(1,1), NH3(2,2), CS(J=l-0), CS(J=2-1) and 
N2H"'"(J=l-0) were applied to estimate gas parameters. For 
brevity, the mathematics of our molecular emission analyses 
is not presented here, but in §[Bl Optical depths were found 
by taking intensity ratios of emission line pairs. Rotational 
(or excitation in the case of N2H^) temperatures could then 
be estimated from LTE assumptions allowing total column 
densities to be calculated. The methods used are only briefly 
described here, but details are presented in references. After 
making some assumptions about core structure, masses and 
densities were estimated. 

Gaussian functions were fitted to all emission lines us- 
ing a minimisation method. Lines with visible hyperfine 
structure (NH3 and N2H^) had each 'satellite' line con- 
strained to be a fixed distance (in VLsn-space) from the 
'main' line. In addition to this, for simplicity, all NH3(1,1) 
emission lines were assumed to have equal FWHM. 



to the centre lines as 'main' lines and all others as 'satellite' 
lines. 

To estimate optical depths, r, of NH3(1,1) and 
N2H'*'(J=l-0), any deviations from the aforementioned rela- 
tive line strengths of main and satellite emission peaks were 
assumed to be due to optical depth e ffects, allowing the us e 
of the standard analysis presented in lBarrett et al. 1 l|l977t ). 
Note that the optical depths quoted for these two molecular 
transitions are those of the central emission lines, which in 
both cases have unresolved structure. 

By assuming abundance ratios between CS and rarer 
isotopologue partners, C'^^S and ^'^CS, a similar method can 
be used to estimate the optical depth of CS emission lines. 
Values of 22.5 and 75 were assumed for the abundance ra- 
tios [CS]/[C-^*S] and [CS]/[^^CS] respectively. CS emission 
without detected isotopologue pairs were approximated as 
being optically thin. 

Where multiple optical depth values were calculated (4 
for NH3(1,1), 2 for N2H+(J=l-0), 3 for CS(J=l-0,2-l) re- 
spectively), an average was taken, weighting by the inverse 
of the optical depth variance. 



4.2 Temperature and Column Density 

Upper-state column densities, Ny, for NH3(1,1), 
N2H+(J=l-0) ar id CS(J=l-0, J=2- l ) we re calculated 
via Equation 9 of IColdsmith fc Langeij 1 19991 ). 



NH3: The LTE an alysis outhned in 



and summarised in lUngerechts et al" 



Ho fc Townel (|l983l ) 
( 19861 ) was used to 



find rotational temperature. Trot, and total column density, 
Niot, of NH3 from the NH3(1,1) andNH3(2,2) emission lines. 
Kinetic temperature, Tkin, was cal c ulated from an approx- 
imation presented in iTafalkT et al. I (|2004l ). 

N2H"'": Given the smaller beam FWHM of the 3mm obser- 
vations with Mopra (36" < 0.2 pc at Ikpc), a filled-beam 
assumption for the measured N2H^(J=l-0) lines was used 
to estimate the excitation temperature, Tex, in order to find 
the total column density as suming an ap proximation of a 
linear, rigid rotator in LTE (|Hotzel et al. | [2004V 



CS: Due to the small seperation in energy (~4.7K), where 
both CS J=l and J=2 column densities were calculated, 
a weighted average was taken and assumed to represent 
the CS(J=1) column density (after a beam-dilution correc- 
tion, see i]4.3p . Where possible, an NH3-derived temperature 
would be assumed, otherwise a temperature of 10 K was as- 
sumed and applied in an LTE approximation to calculate 
the total CS column density. 



4.1 Optical Depth 

The NH3(1,1) inversion transition has 5 distinctly resolvable 
peaks and the N2H"''(J=l-0) rotational transition has 3 (fur- 
ther splitting is present in both molecular transitions, but re- 
main unresolved due to line overlap). The relative strengths 
of these lines (increasing i n velocity-s pace) are [0.22, 0.28, 
1, 0.28, 0.22] for NH3f l.l) JHo fc Tow ncs 198:^) and [0.2, 1, 
0.6] for N2H+(J=l-0) jWomack et al. Il 19921 ). Here, we refer 



4.3 Beam Dilution 

A beam dilution factor, /, and coupling correction factor, K, 
were applied. Multiplying the column density by the factor: 



/i^ = 1 - exp 



4fl" 



In 2 



(1) 



corrects for the size of the so urce compared to the beam 
FWHM l|Urauhart et al. l2010l ). where R is the source radius 
and 6mb is the main beam FWHM. The main beam FWHM 



8 Nigel I. Maxted et al. 



for CS(J=2-1), N2H+(J=l-0), CS(J=l-0) and detected NHg 
lines are 36±3", 36±3", 59.4±2.4" and 12 3±18" respectively 
ijLadd et al. ll2005l : IUrauliart et al. Il201(]h . fK tends to 1 for 
beam sizes smaller than the core size. 



4.4 LTE Mass and Density 

LTE masses, Mlte, were calculated assuming a spheri- 
cal, homogeneous core of radius, R, composed of molec- 
ular hydrogen. Average abundances with respect to H2 
for CS, N2H"'" and NH3 were assumed to be 1x 10"^ 
jFrerking et al. Ill980l). 5><10~ " jPirogov et al. 1 120031 ') and 
2x10"^ jstahler fc Paliall2005l ) respectively. 



4.5 Virial Mass 

Under the assumption that turbulent energy balances grav- 
itational energy within a core, a virial mass: 



R 

1 pc 



1 kms~ 



\Mr. 



(2) 



can be calculated, where R is the radius of the core, Aiifwhm 
is the emission line full-width-half-maximum and A: is a co- 
efficient that depends on the density profile of the core. 
The coefficient k is 4 44 for a Gaussian density profile 
l|Protheroe et al. ||2008| ). 210 for constant density an d 126 
for density, p a R~^ profile (|MacLaren et al.~l Il988l ). Re- 
cen t CO observation s suggest that the latter is the best 
fit (ISano et al. II2OIOI ). so fc=126 is assumed. A core radius 
consistent with previous beam dilution assumptions is as- 
sumed. For CS emission, virial masses were calculated using 
the rarer isotopologue C^*S to minimise the effect of optical 
depth-broadening. 

As the virial mass is not always consistent with LTE 
mass, the molecular abundances that give an LTE mass 
equal to that of the virial mass was calculated in each case. 



5 CORE PARAMETERS 

The critical density of the J=l-0 transition of CS is 
~8x 10* cm~'^, implying that the density of regions traced 
by this emission, in this case. Cores A, C and D, are prob- 
ably distributed around this val ue. Cores A, C and D are 
coincident with infrared sources iMoriguchi et al. Il2005h . so 
the detection of CS(J=l-0) emission from these regions sup- 
ports the notion that these cores harbour dense warm gas. 



5.1 CoreC 

The 3, 7 and 12 mm observations on Core C revealed emis- 
sion from six different molecular species (see Tablel^)). Using 
three of these molecular species, detailed analyses were per- 
formed. Results are displayed in Table |3l and discussed in 
the following subsections. 



5.1.1 CoreC Size 

Using CS(J=l-0) map data, we can estimate the core size 
by deconvolving the intensity distribution with a Gaussian 
of FWHM equal to the 7 mm beam FWHM. We assume 



that the Core C radial density profile is well-approximated 
by a Gaussian and assume the relation = s'^ + r'^, where 
r, s and r' are the observed radius, beam half-width-half- 
maximum and deconvolved radius, respectively. By extend- 
ing the definition of s to be the quadrature sum of both the 
7 mm HWHM (29.7") and the Gaussian smoothing HWFM 
of our maps (37.5 "), we obtained a dense-core radius of 
~25". Thus, we assumed a core radius of 0.12 pc (25" at a 
distance of 1 kpc) in our calculations, yielding correction fac- 
tors, fK, of 1.3, 1.3, 2.5 and 9.1 for CS(J=2-1), N2H+(J=1- 
0), CS(J=l-0) and NH3 lines, respectively. There is a ~50% 
systematic uncertainty associated with our calculated ra- 
dius that propagates into LTE mass (Mlte), density (n^fj), 
virial mass (Alvir) and virial abundance (Xvir). We do not 
account for this systematic in TableO but discuss it in 



5.1.2 Mass and Molecular Abundances 

The adoption of molecular abundance ratios is probably the 
dominant source of systematic uncertainty. It is not uncom- 
mon for molecular abundances to vary by an order of magni- 
tude between different Galactic regions, but given the order- 
of-magnitude agreement between the calculated LTE masses 
(4O-8OM0), the assumed abundances are probably valid. Un- 
der the assumption that turbulent energy balances gravita- 
tional energy within Core C (which may not be valid in a 
core that is shock-disrupted), virial masses (30 to 70 Mq) for 
these three molecular species allowed abundances with re- 
spect to molecular hydrogen to be calculated. We saw no sig- 
nificant deviation in CS, N2H+ and NH3 abundance (~lx, 
~lx and ~2.5x, respectively). 

Systematic uncertainties in beam efficiency, beam- 
dilution and beam-shape corrections axe difficult to account 
for, but propagating published systematic error estimates 
leads to greater than 30% systematic error in N2E[^ and CS 
LTE mass and density estimates, but has only a ~5% effect 
on NHa-derived mass and density estimates, since many such 
effects cancel-out in the analysis. 

In our CS analysis of CoreC, several assumptions were 
made. The use of an NHa-derived temperature introduces 
some uncertainty, as NH3 emission may tra ce a region of 
different temperature to that of CS emission. iTafalla et aTl 
(2004) examined the starless cores L1498 and L1517B in 
detail and found significant variation in CS and NH3 molec- 
ular abundances with radius. NH3 (para-NH3) abundance 
peaked at ~1.4 and 1.7x10^* towards the core centres and 
dropped to ~5xlO"^ at ~0.05pc, while CS (and C^*S) 
abundances decreased towards the centre, with models as- 
suming a central ~0.05 pc hole, likely due to CS-depletion 
onto dust grains, providing a better fit to observations than 
constant-abundance models. It follows that it is feasible that 
the region emitting in CS(J=l-0) and CS(J=2-1) is different 
to that emitting in NH3(1,1). The systematic error caused 
by this effect is difficult to quantify, but given the factor 2 
agreement between CS and NH3-derived core mass under 
a 'perfect-correspondence' assumption, it is expected to be 
small. 

The CS J=l and J=2 states are separated by only 
2.35 K and a comparison of the two emission lines to retrieve 
a rotational temperature would likely be dominated by un- 
certain systematics since the 2 measurements were taken 
with different spectrometers and were separated by 2 years. 
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Table 3. CoreC parameters derived from three different molecular species assuming a core radius 0.12pc. Optical depth, rotational 
temperature, kinetic temperature, molecular column density, H2 column density, LTE mass, LTE density, virial mass and molecular 
abundance are denoted by r. Trot, Tj,i„, Nx, ^LTE, ^H, ^vir and Xvir- Statistical errors (rounded to the nearest significant 

figure) for all variables independent of virial equilibrium assumptions are shown. Systematic errors are not displayed, but are discussed 
in the text. 
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"Equal excitation temperature assumption, ''Assuming a source of radius 25" (~0.12pc), "^LTE assumption, Assuming abundance 
ratios stated in i|4.4l ''Tkin =26±8 ^Tex, ^Assuming the temperature calculated from NH3, ''Using the optically thin isotopologue 

C34S. 



As we expect the J=l and J=2 states of CS to be roughly 
equally-populated, the optical depths and upper-state col- 
umn densities were averaged to minimise instrumental sys- 
tematic uncertainties (deemed to be larger than errors in- 
troduced by the equal-population assumption). 

Of course, all analyses thus far have assumed LTE, 
which may not be valid, as suggested by CS J=l-0 and 
J=2-l intensities. We investigate how the calculated gas pa- 
rameters may differ from a more accurate non-LTE statis- 
tical equilibriu m treatment using the f reely-available soft- 
ware RADEX (|van der Tak et aL II2007I '). Through RADEX 
modeling, we conclude that for a CS column density of 
lxlO^*cm~^, the molecular hydrogen density may not ex- 
ceed ~6xl0''cm~^ at a temperature of 10 K, and may be 
smaller for increasing temperature. This solution would have 
implications for the calculated mass and abundance, be- 
ing ~75% smaller and ~3xl0~^, respectively. We note that 
in order to correctly simulate the observed ^^CS intensity 
using the gas parameters outlined above (which were con- 
sistent with CS and C'^^S emission), the assumed value of 
[CS]/["CS] must be adjusted to ~110. Although a limited 
investigation of CS isotopologue abundance is presented in 
i^ lA41 we leave a more detailed examination of this for later 
work. 

We note that our measurements see m to trace smaller 
mass than CO(J=1-0) observations by iMoriguchi et al.1 
l|2005ll , who calculate a Core C mass of ~400 M© . This sug- 
gests that CS(J=l-0), NH3 and N2H^ are tracing different 
regions of the molecular cloud to CO, most likely a denser, 
inner component. The Core C radial density profile is dis- 
cussed further in §|A3l 



5.2 Cores A and D 

Table [4] shows gas parameters derived for Cores A and D. 
We estimate the Core A radius to be ~0.3pc (60" at Ikpc) 
via the method described in i]5.1.1l Core A, being more ex- 
tended than Core C, and having an infrared source possibly 
offset from the CS(l-O) peak emission (see i) l3.1|) . may be 
composed of two smaller clumps including one emitting in 
infrared. These imhomogeneities may even be the result of 
an outfiow originating from IRAS 17082—3955, but this sce- 
nario remains as speculation. As our maps do not encompass 



CoreD, we cannot estimate its radius and assume a radius 
of 0.12 pc for CoreD (see ij |A2l for alternate radius assump- 
tions). 

Given our assumptions, LTE and viral mass solutions 
for Cores A and D range from 12 to 120 Mq . We calcu- 
late CS abundances a factor 0.2 to 1 of that assumed, 
and, like in CoreC, we note that CS(J=l-0) seems to 
trace less of the total mass than CO(J=1-0). Core A and 
D CO-derived masses are ~700 and ~3OOM0 respectively 
(Moriguchi et al. 2005). Similar to CoreC, CS likely traces 
the denser, inner components of Core's A and D. 



6 HIGH ENERGY PARTICLE PROPAGATION 
INTO DENSE CORES 

In the context of understanding the nature of the TeV 
gamma-ray emission from SNRs, of considerable interest is 
the energy-dependent propagation of high energy protons 
(CRs) and electrons into dense cores, which may lead to 
characteristic gamma-ray spectra from GeV to TeV ener- 
gies. This effect could offer a new way to identify the parent 
particles (hadronic vs. leptonic) responsible for gamma-ray 
emission, pr ovided such cores ca n be spatially resolved in 
gamma-rays. Gabi ci et al. ](l2007') investigated the propaga- 
tion of CR protons into molecular clouds and found that 
in some models, suppressed diffusion of low energy CRs re- 
sulted in a gamma-ray spectrum towards the cloud centre 
that is harder than that towards the cloud edges. This effect 
became more pronounced for more centrally-condensed gas 
distributions. Similar phenomena may in fact be occurring 
towards the dense cores associated with RX J1713.7— 3946. 

We therefore investigated CR diffusion into a molecu- 
lar core located adjacent to a CR accelerator. Our model 
is a simplistic approach that only aims to show in principle 
the effects that dense gas could have on the spatial distri- 
bution of CR protons of different energies, which will have 
knock-on effects on the gamma-ray and X-ray spectra from 
various parts of the core. The model configuration is based 
on the Core C molecular core, and its apparent location with 
respect to RX J1713.7— 3946 as an assumed CR accelerator. 
CR transport in this region is expected to be in the diffusion 
regime, since CR (proton) gyroradii at TeV energies are ex- 
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Table 4. CS(J=l-0) parameters for Cores A and D assuming a kinetic temperature of 10 K. No emission from a rarer isotopologue of 
CS is detected for these cores, so the maximum optical depth is calculated from the ITrms noise in the C^^S band. We display derived 
values for both optically thin and thiclc assumptions for parameters that depend on optical depth, giving a range of solutions. Optical 
depth, molecular column density, H2 column density, LTE mass, LTE density, virial mass and molecular abundance are denoted by r, 

Ncs. Nh2. Mlte, i^h, Myir and Xvir- 



Core 


T 


Ncs 




Mlte 






'X.vir 






(xl0i2cm-2) 


(xl02icm-2) 


(Mq) 


(xlO^ cm-3) 


(Mo) 


[CS]/[H2] 


Core A 


0-0.44 


3-4 


3-4 


12-15 


3-4 


60 


2-3x10-1'^ 


CoreD 


0-4.1 


40-170 


40-170 


30-120 


100-500 


100 


3-10x10-10 



"■Assumed rotational temperature of 10 K, ^Corrected for beam dilution assuming a source of radius of 60" (^0.30 pc) for Core A and 

25" (^0. 12 pc) for CoreD, ''Assuming abundance ratios stated in ^4-4\ 



pected to be <10~*pc for typical magnetic field strengths 
of 10 fiG or more. 

Our model core is placed 5pc from the centre of 
RX J1713.7— 3946, which has an outer shock radius growing 
with time to eventually reach the core, as in the situation 
assumed for Core C. We then calculated the propagation dis- 
tance (via diffusion) of CR protons of energy, Ep, from the 
outer shock of the SNR, towards the centre of this molecular 
core, applying time limits based on assumptions concerning 
the age of RX J1713.7— 3946, and the escape time of CRs 
from the SNR outer shock. 

The Molecular Core: We modelled a core with radius set 
equal to 0.62 pc and a density of 300 cm~^, consistent with 
volume-avera ged density estimates from C0(J=2-1) mea- 
surements by ISano et al. for Core C. The diffusion 
of protons is dependent on the strength of magnetic turbu- 
lence within the region. These magnetic fields are generally 
assumed to be frozen-in to the gas, such that denser gas im- 
plies larger magnetic field strengths, -B(nH,). If we assume 
the relation stated in ICrutcheil i|l999l ): 

i?(ng,) ~ 100/ "^^3 [mG] (3) 
V 10 cm ^ 

the average magnetic field strength inside the core is 

B r^l7flG. 

6.1 Cosmic-Ray Proton Transport 

In the hadronic scenario for gamma ray production, CR pro- 
tons diffuse from acceleration sites and collide with matter 
to create TeV emission. The rate of diffusion of CRs is gov- 
erned by magnetic turbulence, the characteristics of which 
are somewhat uncertain. 

If we assume t hat RX J1713.7-3946 is ~1600 years old 
l|Wang et al. Ill997l '). we have a limit on the amount of time 
CR protons have had to be accelerated in the SNR shock 
and diffuse into the core. Our goal was to examine how far 
CR protons of different energies penetrated into the core. 
In Figure fTUl we have visualised the results (a slice view of 
the core) of the following equation, which gives the aver- 
age radial penetration distance of a CR proton within our 
modelled molecular core. 

R = 0.62 - ^6D{Ep,B)[1600 -to] [pc] (4) 

where D{Ep, B) is the diffusion coefficient that depends on 
CR energy, Ep, and magnetic field strength, B. Equation[4] 



returns the average distance a CR of energy, Ep, will pene- 
trate into the modelled molecular core given 1600 years. 

The start time, to, was set equal to the minimum of the 
average time taken for the CR proton to escape the SNR 
shock, tesc, and the time taken for the SNR shock to reach 
the boundary of the core fcore- 

to = min(fesc, tcore) [yr] (5) 

The value for tcore ~600 years was based on the Sedov so- 
lution for a core-collapse supernova explosion into a wind- 
driven bubble, w here the shock radius evolu tion is pro- 
portional to t"''/* (|Ptuskin fc Zirakashviiill2005l '). The time- 
dependent shock radius and energy-dependent escape time 
mean that some CRs will escape the shock before it contacts 
the core. After the shock has passed the core (~600 years), 
it is assumed that CRs of all energies can reach the core 
boundary. In general however for the moment we ignore this 
CR diffusion time from the shock to the core, which is at 
most ~100yr for the highest energy CRs since it is assumed 
to take place in a relatively low magnetic {B ~few /iG) and 
density (n < Icm"^) environment. 

Escape Time: With the assumption that within a SNR 
the highest energy particles (3>lTeV) are accelerated in 
the tr ansition from the free-expan sion phase to the Sedov 
phase (^Pt uskin fc Zi rakashvili' 2005) within a short timescale 
( Uchivam a et al. 112007), wc used the f ollowing parameter i- 
sation for the CR proton escape time, (|Gabici et al. |[2009l l: 

tesc{Ep) = tsedov { I [yr] (6) 

\liP,max J 

where Ep^rnax is the maximum possible CR pro t on en ergy, 
set equal to 5 x lO'^^ eV following ICaprioh et aTl (|2009l ) and 
'Cas anova et al. I(|2010l ). tsedov is the time for the onset of the 
Sedov phase, set equal to 100 years, and 5 is a phenomono- 
logical parameter that describes the energy-dependent re- 
lease of CRs. 5 is set equal to 2.48 assuming that ~10^^ eV 
CRs are accelerated near the start of the Sedov phase re- 
ducing to ~10® eV near the end (see Gabici et al. 2009). 

Diffusion: ICabici et alTI l|2007h considered the timescales 
involved in different CR transport processes within a 
dense 20 pc-diameter molecular cloud. They demonstrated 
that the dynamical (free-fall) and advective (turbulent) 
timescales involved were likely an order of magnitude larger 
than the time-scale for inelastic proton-proton interactions 
and a varying amount larger than diffusion time-scales. 
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Figure 10. The logarithm of the minimum energy CR proton 
able to penetrate into different radii of our simulated core. The 
results, presented as 2D slices, from 3 different diffusion suppres- 
sion coefficients, x ^^e displayed. Note that the arcsec scale for 
penetration distance assumes a distance of 1 kpc. Hatching indi- 
cates a lack of CR inhabitanco of CRs. 



For our model core, we calculated dynamical and advective 
timescales of >10* years, so we too ignore these mechanisms 
and only consider proton-proton interactions (see later) and 
difltusion, using a difltusion coefficient parameterised as: 



D{Ep,B{r)) = xDo 



Ep /GeV y 
B/3fiG J 



(7) 



where Do is the galactic diffusion coeffi cient, assumed to 
be 3 x lO^'^ cm^ s~^ to fit CR observations l|Berezinskii et alTI 
Il990l l. and X is the diffusion suppression coefficient (assumed 
to be <1 inside the core, 1 outside), a parameter invoked to 
account for possible deviations of the av erage galactic dif- 
fusion coefficient inside m olecular clouds (jBerezinskii et alTI 
Il990l :[ Gabici et al. |[2007f ) which is largely unknown. 



Interactions: For energies above ~300 MeV proton-proton 
p-p) interactions dominate over ionisation interactions 
Gabici et al. |[2007l ). so we consider the rate of hadronic CR 
interactions with molecular gas protons and ignore ionisa- 
tion losses. Assuming an inelasticity of 0.45 (2 interactions 
give ~79% energy loss) and a proton-proton cross section 
of 40 mb, appropriate for TeV energies, the lifetime of this 
proton-proton interaction process is; 



X 10 



5 f nnir) 



\ 100 ( 



(8) 



Thus, CRs in the core have a large proton-proton interaction 
lifetime of ~10^ yr, so we do not model this process. 

Diffusion Suppression and Results: Diffusion- 
suppression coefficients x = 0.1, 0.01, 10"^, 10"" and 
10~^ were trialled for proton energies 10^^°, 10^"'^, up 
to 10" eV. This wide range of x values were chosen to 
investigate the level of CR penetration with results from 
the smaller x values displayed in Figure fTOl 

For x=0-l ^nd 0.01, all CR energies were able (on av- 
erage) to reach the inner core within the age of the system. 
For x=10~^, some energy-dependent penetration could be 
seen, as CRs of energies 10^" and 10^"'^ eV did not reach the 
inner core. Energy separation became more prominent with 
smaller values of x, as illustrated in Figure fTOl For x=10~" 
and x=10~^, CRs of energy 10^^'^ to > 10^** are now ex- 
cluded from the central regions of the core. This effect would 
have implications for the spectrum of gamma-ray emission, 
with the TeV gamma-ray spectrum becoming progressively 
harder towards the core's centre. 

It is unclear whether such extreme cases of diffusion- 
suppression down to x=10~^ required to produce 
these energy-dependent penetration effects on such small 
scales are plaus i ble. R adio continuum measurements by 
iProtheroe et al.1 l|2008l ) when compared to the GeV gamma- 
ray emission have suggested upper limits for x < 0.02 to <0.1 
(dependent of the assume d magnetic fie l d) ins ide the Sgr B2 
molecular cloud (see also I Jones et al.1 l|201ll )V Suppressed 
CR diffusion (x ~ 0.01 to 0.1) is also indicated when ex- 
plaining the GeV to TeV gamma-ray spectral variation seen 
in the vicinity of several evolved SNRs (e.g. 'Gabic i et al. I 
(|2010t )). Moreover it has been recognised that suppressed dif- 
fusion could be expected in turbulent magnetic fields likely 
to be found i n and around shock-disrupted regions (e.g. 
lOrmes et al.1 (|l988D ). Overall, we regard these uncertain- 
ties in particle transport properties as motivation for inves- 
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tigating the efTects of a wide range of diffusion suppression 
coefficients, as applied to our model core which may also be 
similarly shocked and turbulent. 

From a more theoretical standpoint, the level of CR 
diffusio n suppression that may be expected is still somewhat 
unclear. lYan et al. I (l20f 2h recently investigated the effect of 
CR-induced streaming instabilities, distortions of magnetic 
field lines caused by the flow of CRs, on CR acceleration 
within SNRs. The authors expressed interest in solving for 
a diffusion coefficient in the region surrounding SNRs by 
incorporating both streaming instabilities and background 
turbulence. Their future work will help to constrain the level 
of CR diffusion-suppression that can be plausibly expected 
inside gas associated with SNRs. 

We also note that an additional CR component from the 
Galactic diffuse background 'sea' of CRs will be incident on 
the core. Given this is an ever-present component, the pen- 
etration of these diffuse CRs will not be limited by SNR 
age, but by the age of the core, and possibly by the energy 
loss timescale r^p due to proton-proton collisions (Eq|8]). 
For our model core with average density n =300 cm"'', 
Tpp will not be the dominant loss mechanism, but will be- 
come more significant in the central region (r^p ~600 yr for 
n =10^ c m~^). The dynamical a ge of Core C is estimated at 
~f 0^ yr (jMoriguchi et al."1l2005l ) which is similar to the en- 
ergy loss timescale for our averaged-density model core. In 
most of the x values we consider, these timescales are signif- 
icantly larger than the time taken for CRs to diffuse into the 
core centre. For an example extreme case, it takes ~5x 10^ yr 
for a 10^° eV CR to reach the core centre for x =10~^, so 
diffuse CRs will likely penetrate the core. However, the en- 
ergy density of diffuse CRs, assumed to be approximately 
similar to the Earth-like level of CRs, would be consider- 
ably smaller (a factor ~fO~*) than the local CR component 
due to an adjacent SNR (e.g. Aharonian & Atoyan f996), 
so we can therefore neglect it here. 

Our core model reflects only average values of density 
and magnetic field as so will have limitations in comparison 
to a core where these quantities vary with radius. Compared 
to the case for a varying diffusion coefficient resulting from a 
radially dependent density (see Equations [3] and [7]) , our ho- 
mogenous core model may under-estimate the CR penetra- 
tion depth in the outer region of the core, and over-estimate 
the penetration depth towards the inner region of the core. 
In future work we will consider CR propagation into a core 
with a radially varying diffusion coefficient based on the 
measured power-law radial density profile fro m CO mea- 
surem ents n(r) ~ 10^^/(1 -I- (r/O.fpc)^'^) cm"^ ()Sano et al. I 
I2OIOI ). This model should be three-dimensional and predict 
a line-of-sight view of gamma/X-ray emission morphology 
towards the core. 

Another limitation concerns the implications of the ex- 
tremely low values of diffusion suppression we have used 
down to X = 10~^. For such cases, which represent highly 
turbulent magnetic fields, the possibility of second order 
Fe rmi reacceleration of th e CRs might apply as d iscussed 
bv lDogiel Sharovl (jf99Gl l . lUchivama et al. I (|2010l l have in 
fact considered this effect in accounting for the GeV to TeV 
gamma-ray spectra towards several evolved SNRs (W5fC, 
W44 and IC443) by considering the reacceleration of dif- 
fuse CRs inside clouds disturbed by a supernova shock. 
Such effects may occur in molecular cores associated with 



RXJ1713. 7-3946, particularly those contacted by the SNR 
shock but at this stage we leave this as an open question. 
Finally, there may be a regular (or ordered) magnetic field 
component, which if present may either assist or hamper the 
penetration of CRs depending on the strength and orienta- 
tion of this field component. Such modeling is beyond the 
scope of this investigation, but should be considered in more 
detailed models. 



6.2 Gamma-rays from Core C 

We can expect broad-band GeV to TeV gamma-ray emission 
from CRs penetrating the core as they interact with core 
protons. The level of the gamma-ray emission will reflect 
the energy and spatial distribution of CRs within the core 
coupled to the mass of protons they interact with (see for e.g. 
iGabici et al.1 (|20o3)). Variations in the gamma-ray spec- 
trum will reflect the energy dependent penetration depth of 
CRs, and would be a feature only resolvable in gamma-rays 
with an angular resolution of f arcmin or better according 
to the size of the cores. The ~6 arcmin angular resolution 
offered by If.E.S.S. is at present not sufficient to probe for 
such spectral variations. 

For a wide range of diffusion suppression factors x = 
fO"^ to fO"^ our simple model of Core C suggests that lower 
energy CRs can be prevented from reaching its centre. For 
the more severe values oi x = 10~* and fO~^, the CR pen- 
etration to the core centre can be limited to CRs of TeV or 
greater energies, which could lead to considerable differences 
in the gamma-ray spectra seen towards the centre and edge 
of the core. 

The overall gamma-ray fiux expected from dense 
cores towards RX Jf7f 3 .7 — 394 6 has been discussed by 
IZirakashvih fc AharonianI (j2010h (see their Fig.f4). They 
suggested an enhanced TeV hadronic component towards 
molecular cloud cores (with combined mass ~f 000 Mq) may 
be found at a level of ~6x fO~^^ ergcm'^^s"^ for energies 
above ~fOTeV, which would result from CRs of energy 
above ~fOOTeV entering the core, which could be the case 
for most levels of diffusion suppression. For these gamma- 
ray energies the dense-core hadronic component may reach 
or even exceed the SNR-wide leptonic inverse-Compton com- 
ponent. The >fOTeV gamma-ray flux from the central re- 
gions of Core C, which is 40 to 80 Mq based on our measure- 
ments, may therefore reach a fewx fO"^'' ergcm^^s^^, pos- 
sibly making it detectable by future more sensitive gamma- 
ray telescopes such as the Ch erenkov Telescope Array (CTA) 
l)The CTA Consortiumllioiol l. With its expected 1 arcmin or 
better angular resolution, CTA could also potentially probe 
for energy dependent penetration of the parent CRs, by ex- 
amining the flux and spectra from the inner to outer regions 
of the core. 



6.3 Electron Transport and X-ray Emission 
Towards Core C 

As discussed in ijl] two synchrotron X-ray emission peaks 
can be observed towards the boundary of Core C (Figure 
[TJj). These peaks are perha ps due to increased s hock interac- 
tion with the dense gas (see lSano etHI (j2010h :ll noue et al. I 
(20f2)) giving rise to a region of enhanced magnetic field 
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and/or electron acceleration. The observed synchrotron X- 
ray emission can be produced either by electrons accelerated 
directly by the SNR shocks outside of Core C or, by the sec- 
ondary electrons from CR proton interactions with ambient 
gas. In the latter case, the X-ray emission spectrum could be 
influenced by the energy dependent penetration of CRs into 
the core and for strong CR diffusion suppression (x ^ 10"'*), 
one might expect an X-ray spectral hardening towards the 
core centre, in the same way as for gamma-rays discussed in 

The dense gas and high column density from Core C will 
also suppress the X-ray flux for energies less than ~5 keV via 
photoelectric absorption. Thus some X-ray emission may ac- 
tually originate from within or behind CoreC. The X-ray 
flux decrease seen in going from the boundary of Core C to 
its centre is ~95%. Assuming the cross section for photoelec- 
tric absorption by hydrogen, ~10~'^''cm'^ (at ~lkeV), 
a flux decrease of 95% corresponds to a traversed column 
density of N^f ~3x 10^^ cm~'^, similar to the column den- 
sity towards CoreC from our CS, N2H''' and NH3 observa- 
tions which pertain to the core centre (See Table[3]). There- 
fore it is possible that keV X-ray emission, either from sec- 
ondary electrons produced deep within Core C from pene- 
trating CRs or from external electrons directly accelerated 
outside and then entering Core C, are absorbed by gas within 
the core. 

A potential way to discriminate between these two sce- 
narios would be to have arcmin angular resolution imaging of 
>10keV X-rays, for which photoelectric absorption is neg- 
ligi ble. Forthcoming hard X-ray telescopes such as Astro - 
H dTakahas hi et al. I I2OI0I ) and NuSTAR (|Harrisonl I2OI0I ') 
are expected to have such performance. Any >10keV X- 
ray emission centred on or peaking towards the centre of 
Core C may indicate the presence of penetrating high en- 
ergy CRs, particularly if the energy spectrum of the X-rays 
hardens towards the core. If this effect results from sup- 
pression of CR diffusion as discussed in i; l6.f I the same sup- 
pression effect will apply to external electrons accelerated 
outside the core. Strong synchrotron radiative losses on the 
increasing magnetic field towards the core centre may fur- 
ther limit the X-ray component and/or alter the spectrum 
from these external electrons, to the effect that X-ray emis- 
sion may appear somewhat edge- or limb-brightened in con- 
trast to a more centrally peaked component from secondary 
electrons inside the core. For example the synchrotron cool- 
ing time tsync ~ 1.5(B/mG)~*'^ {e/keV)~°-^ yr of electrons 
of energy E responsible for X-ray photons of energy e where 
E « 62.5Y/(e/keV)/(B/mG)TeV, can be compared with the 
diffusion time td ^ 16x~*(-R/lpc) V(B/3MG)/(E/GeV) yr 
into the core with radius, R — 0.62 pc. For hard X-ray pho- 
tons e ~ lOkeV and suppressed diffusion x < 10~^ , we find 
that taync Can become similar to or considerably less than 
td as B increases beyond > 80/iG towards the centre of the 
core. For smaller values of x, isync becomes less than td for 
even smaller values of B. 



7 SUMMARY AND CONCLUSION 

We used the Mopra 22 m telescope to map in 7 mm lines a 
region of gas towards the western rim of the gamma-ray- 
emitting SNR RX J1713.7— 3946, encompassing the molecu- 



lar cores A, B, and C (labeled bv lMoriguchi et al. I l|2005l )'). 
Deep 7 mm and 12 mm pointed observations were also taken 
towards several of these regions (including Core D in the 
eastern rim) and we analysed archival 3 mm data towards 
Core C. Our goals were to investigate the extent and proper- 
ties of the dense gas components towards RX J1713.7— 3946 
and to complement previous CO studies of the low to 
modcratcly-dcnsc gas_ jFukui et al. II2003I : iMoriguchi et al. I 
l 2005. : ,Fukui 2008 : Saiio et al. Il2010t ). The major conclusions 
from our observations are summarised as follows: 

1) Detection of CS(J=l-0) emission towards three of the 
four cores sampled (Cores A, C, and D which have coincident 
infrared emission) confirming the presence of high density 
(> Wcm-^) gas. 

2) Detection of moderately broad CS(J=l-0) emission 
towards a position in between Cores A, B and C (labeled 
Point ABm) which is towards the X-ray outer shock. This 
broad gas may result from the SNR shock passing through, 
although no shock-tracing SiO emission was observed at this 
position. 

3) The LTE mass for the brightest core, CoreC, was 
estimated using three different molecular species, with re- 
sults ranging from ~40 Mq from CS observations to 80 Mq 
from NH3 and N2H"'" observations. The range of masses is 
most likely attributed to uncertainty in molecular abun- 
dances with respect to molecular hydrogen. Virial assump- 
tions allowed molecular abundances with respect to molec- 
ular hydrogen of ~5xl0~*, ~5xl0~*° and ~0.7 to 1x10"'' 
for NII3, N2H"'" and CS respectively to be calculated for 
Core C. Preliminary non-LTE RADEX modeling suggests a 
factor 2 lower density than our LTE results 

4) Cores A and D were found to have LTE masses of 
12 to 15 M0 and 30 to 120 M0 , respectively, as traced by 
CS(J=l-0) emission, although these values are subject to 
core-radius uncertainties. 

5) The inferred column density of our CS and N2H"'" 
measurements towards Core C could lead to considerable 
photoelectric absorption (~95%) of the < 5 keV X-ray emis- 
sion. The small-scale keV X-ray features on the border of 
Core C therefore might not represent the complete X-ray 
emission towards this core. 

6) We investigated energy-dependent diffusion of TeV 
cosmic-ray protons into Core C using a simple model for the 
core based on average values of density and magnetic field. 
For the cases of suppressed diffusion coefficients with factors 
~ X =10"^ down to 10~^, lower than the galactic average, 
considerable differences in the penetration depths between 
GeV and TeV CRs were found, with GeV CRs preferentially 
excluded. This effect could lead to a characteristic harden- 
ing of the TeV gamma-ray and hard X-ray emission spec- 
trum peaking towards the core centre. Such features may 
be detectable and resolvable with future gamma-ray and 
X-ray telescopes and therefore offer a novel way to probe 
the level of accelerated CRs from RX J1713.7— 3946 and the 
strength and structure of the magnetic field within the core. 
A more detailed diffusion model assuming a variable density 
and magnetic field is left for later work. We would also finally 
note that mapping i n tracers of cosmic-ray ionisation such as 
DCO+ and HCO+ (|Montmerlell2010l : llndriolo et al. Il2o"loh 
could a be highly valuable link between RX J1713.7— 3946 
and the observed gas and provide complementary informa- 
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tion on the level of low energy MeV to GeV cosmic-rays 
impacting the molecular cores. 
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APPENDIX A: THE MOLECULAR CORES IN 
FURTHER DETAIL 

Al C^^S and HC3N emission 

C^*S(J=l-0) emission can be viewed towards CoreC (Fig- 
ure lAll top). This is further evidence of the existence of a 
dense molecular core and, with the CS(J=l-0) spectral line, 
allows the calculation of the optical depth and column den- 
sity. HC3N(J=5-4) emission from Cores A and C (Figure|XTl 
middle and bottom) are also indicative of dense cores. 



A2 Scaling Factors 

In this investigation, the CoreC radius of 0.12pc, as cal- 
culated from CS(J=l-0) emission, has an uncertainty of 
~50% and may not necessarily be representative of core sizes 
traced by NH3(1,1) and N2H"'"(J=l-0) emission. This radius 
was also assumed for CoreD. We therefore provide scaling- 
factors (Table RT|) normalised to a core radius of 0.12 pc to 
allow parameter estimates derived from the three molecular 
species to be readily scaled to represent different source size 
assumptions. These were derived by simply calculating gas 
parameters for other core radius assumptions and dividing 
this by the result for a core radius of 0.12 pc. 



A3 Core C Radial Density Profile 

ISano et al. I l|2010h used CO measurements to derive a 
power-law radial density profile of index — 2.2±0.4 (Figure 




Galactic Longitude (deg) 

Figure Al. Top: C3-*S(J=l-0) vlsr = -13.5 to -9.5 km s"!. 
Middle: HC3N(J=5-4) vlsh. = -13.5 to -9.5kms-i. 
Bottom: HC3N(J=5-4) vlsr = -11.0 to -9.0kms-i. 
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Table Al. Mass, M, density, n^f, and abundance, Xvir, scaling 
factors versus assumed CoreC radius. Gas parameters derived 
from an assumed radius of 0.12 pc can be readily converted to 
be for other assumed core radii by multiplying the value by the 
appropriate correction. 



Core Radius (pc) 


0.05 0.1 0.15 0.2 


0.25 


0.3 


M 


Mass Scaling Factors 






LTE 12 mm 


5.7 1.5 0.68 0.41 


0.28 


0.22 


LTE 7 mm 


5.3 1.4 0.72 0.50 


0.41 


0.37 


LTE 3 mm 


3.7 1.2 0.85 0.76 


0.75 


0.74 


Virial 


0.42 0.83 1.25 1.67 


2.08 


2.50 




Density Scaling Factors 






LTE 12 mm 


83 2.7 0.37 0.09 


0.03 


0.01 


LTE 7 mm 


77 2.6 0.39 0.11 
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0.03 


LTE 3 mm 


53 2.3 0.46 0.17 
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\vir 


Molecular Abundance Scalin 


g Factors 


12 mm (NHs) 


14 1.8 0.56 0.25 
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13 1.7 0.59 0.30 
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Figure A2. Ave rage density within radius, as traced by CO 
dSano et al. II2OI0I) and CS (this work). A —2.2 power-law func- 
tion is also displayed for comparison (not fitted). Uncertainties 
are not displayed and are discussed in the text. 



Extrapolating the model to R <0.12pc yields a num- 
ber density value of ~2x 10* cm"''. Our LTE results indicate 
a density of 2 to SxlO^cm""^ and virial densities are 1 to 
4x10* cm~'^, all an order of magnitude larger than that in- 
dicated by CO. 

It follows that the Core C density profile may increase 
beyond the CO-derived power law at smaller radii, but this 
discrepancy could be due to abundance uncertainties dis- 
cussed in H5.1.2I This would not account for the discrep- 
ancy in virial-derived densities, which may simply be due 
to the virial assumption being invalid for Core C. If the ap- 
parent core radius varies between different molecular tran- 
sitions, densities derived from different emission lines could 
be scaled independently (see Table f^TI) . potentially changing 
estimates by an order of magnitude. 



A4 CS Isotopologues 



A key assumption 
lations is that of 
tope ratio values, 

rl3 



in the optical depth calcu- 
the sulphur and carbon iso- 
[CS]/[C=*''S]~p2S]/p''S]~22.5 and 
[CS]/['^CS]~[i2C]/['^C]~75, respectively. These assump- 
tions can be investigated by comparing the emission of 



a The Sickle (GO HII O Emission Nebula 
region); Lower lim- 
its 

□ Other SNRs 

■ RXJ1713, Core C 
{This work) 

□ Dust ridge 

■ Open Clusters & 
Star-forming 
Molecular clouds 



□ HII region 

M starburst galaxy 

□ GO (SgrA and B) 







0.05 0.25 0.45 0.65 0.85 1.05 1.25 1.45 1.65 1.85 2.05 2.25 



^^CS/C''''S integrated intensity ratios for the 
and J=3-2 transitions. Data taken from 



Figure A3. 

J=l-0, J=2-l 

iFrerking et al. 1 lll980h. iMartin et al. ] liooi), I Martin et al. I 
( 2005''l and lNicholas et al. Il|2011^ 7 Typical statistical error in data 
values is 25%. 



"CS and C^^'S. The "CS(J=l-0)/C='''S(J=l-0) peak 
intensity ratio from Core C is 0.26±0.06, which would 
appear to be consistent with the adopted abundance 
ratios (22.5/75=0.3). On the other hand, the integrated 
intensity ratio is 0.16±0.05, somewhat inconsistent with 
our abundance assumptions. This is of concern as inte- 
grated intensity is more relevant than peak intensity when 
calculating molecular abundances. We reviewed ^'^CS(J=1- 
0)/C^*S(J=l-0) integrated intensity ratios in published 
literature towards a variety of cores. 

Figure lXSl presents a histogram of the ratio of integrated 
intensities of ^^CS and C^^S emission from this sample. The 
graph suffers from low number statistics, but it can be seen 
that there is a large range of recorded ["CS]/[C^''S] val- 
ues, and although the CoreC RX J1713.7— 3946 data point 
is at the low end of the distribution, it represents only a 

deviation from the average value of 0.55 (excluding the 
Sickle sources). Note- worthy in Figure IXSl are two galactic 
centre sources that recorded [^''CS]/[C^''S] values of over 4(t 
from the mean. More data would be useful to compare the 
["CS]/[C^''S] ratio of CoreC to that of other cores. 



APPENDIX B: THE LTE ANALYSES IS DETAIL 
Bl Ammonia Analysis 

Bl.l Optical Depth 
TX{J,K,m) _ 1 



-T(,/,X,m) 



TX{J,K,s) 1 
t{J, K, m) 



tot 
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B1.2 Rotational & Kinetic Temperature 



Trot = -41 



In 



9 T2 o 



20 tJ°i* 



[K] 
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B1.3 Column Density 
91,1 

Q{T) = Siffie-^*/'^^ 



tot \ p 

Tmbdv [cm-2](B6) 

1 - e ^1.1 / J 



^^E,,,/kTQ^j,^ [cm-^] 
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B4 LTE mass, Density and Virial Abundance 



M = 



2-Kr tuhNmoi 



Xvir — 
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Xr 

2Trr'^mHNMoi 
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B2 N2H+ Analysis 

B2.1 Optical Depth 



Ttot = 



1 - e- 



(B9) 
(BIO) 



B2.2 Excitation Temperature 

m Tmb . I AOs 



Tb 



(Bll) 



Column Density 
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B3 CS 

Optica/ DepJ/i 
Tcs 1 - e-^cs 



Tcs zso. 1 - e-"^cs 



(B13) 



where a is the ratio of the rarer CS isotopologue abundance 
and CS abundance. 



B3.2 Column density 

^^^^^^cs(j=i)_^Ecs(j=^)/kTQ^j.^^^^ [cm"'] (B15) 

gcs(j=i) 



B3.3 Rotational Temperature 
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